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INVESTIGATION OF RESONANCE INTEGRALS OCCURRING 


IN COSMIC RAY DIFFUSION THEORY 


We have investigated the method of calculating the 
pitch angle diffusion coefficient for cosmic-rays in a static 
random magnetic field using the "resonance integral" method 
of Hasselman and Wibberenz [1968]" and Jokipii [1972] . As 
elaborated by Jokipii [1972] the pitch angle diffusion 
coefficient may be derived from the Vlasov equation via 
ensemble averaging to obtain the result 


D(y) 


dkP (k) 


sin (k-k Q ) z 
k-k 


(1) 


where k Q = (jJ 0 /yw, z = ywt,u) 0 is the gyrofrequency in the 

average field <B>, u is the cosine of the pitch angle and w 

is the particle speed. P (k) is the wave number power spectrum 

of the random part of the magnetic field. 

Jokipii presents arguments that for t>>o^ 1 the resonant 
sin (k-k ) z 

function — becomes effectively a delta function of 

k-k 

o 

k and may be replaced by TrS(k-k 0 ) to give 


W 

D(y) ^(t^) 

MW 


(2) 


We have investigated the integral in equation (1) by means of 


contour integration using a power spectrum of the form 


P(k)=k~ [!+(£-) 2 ] 01 where k = % ~ l , o 


being the correlation 


- i 
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length, and j <a<l. We obtain the integral 

2 a-i r°° i(k-k Q )z -i(k-k Q )z 

I = -£ / dk § — (3) 

2i J (k+ik c ) a (k-ik c ) a (k-k o ) 

Since there is no pole at k = k Q for the total integral 

we may run the contour above pr below the real k axis provided 

we do so systematically. After this choice has been made 

the integral may be split in two and the contour for the 

integral containing e o' (e o ) may be moved up, 

c (down, c ) deforming about the (now present) pole at 
1 2 

k = k Q and the branch points at k = ±ik c as in figure 1, where 

it is assumed that the initial contour ran below k = k . 

o 

The pole term yields the time independent tern 

- Jr ti + (e 2 ) 2 !' 0 - **<k 0 > < 4 > 

c c 

which is the "resonance" result of Jokipii and Hasselman and 
Wibberenz. The integrals along the branch cuts may be trans- 
formed into real integrals and numerically integrated. The 
result is given in terms of the ratio of time dependent to 
time independent terms R(e,i) where the total integral is 
given by I = I Q [1 + R(e,f)] and where e = k c /k Q and 

t = k Z = Z/ & . Our result is 
c c 
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R(eft) = — sin-rra e a (l+e 2 ) a e T 

TT 


1 - 2 , 


x [sin( T /e)F (e , t) - ecos (x/e) H (e ,x ) ] 


(c,I)= $ 


dx 


»-xx 


(x 2 +2x) u [l+e z (1+x) ~] 


( 5 ) 


H (e fx) = F(e,x) - 


F(e,x) is a slowly varying function of x of order unity and 

. - 1 -2 

is plotted for e= 10 , 10 and a = 3/4 in figure 2. 

We see from equation (5) that the ratio is large 
(e 1 2a >>l for a> j) / oscillating at the gyrof requency 
(x/e = w Q t) , and decays exponentially on the correlation time 


scale x = Z/£ not in a time to 
c o 


-1 


It would appear that for 
this form of the diffusion coefficient the arguments of Hassel- 
man and Wibberenz and of Jokipii do not hold. 

These authors also calculate the ad hoc Fokker Planck 
coefficient <Au 2 >/At by calculating <[Ay(t)] 2 > as a function 
of time, dividing by the time t and demonstrating that this 
ratio approaches a constant asymptotically as t-* 00 . This 
constant is equal to the asymptotic value of D(u) obtained 
from the Vlasov equation. The expression for < [Ay (t) ] 2 > is 


< [Au (t) ] 2 > <*/ dk 


p (k) [1-cos (k-k Q ) z] 
__ _ 


<k-k o )* 


( 6 ) 


where all of the symbols have the same meanings as before. 

The resonance argument now leads to the replacement of 
the term [1-cos (k-k Q ) z] / (k-k Q ) 2 by uz6 (k-k ) and yields 


< [Ay (t) ] 2 > oc 7 ttP(k ) 

o 


( 7 ) 
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Since the remainder terms in D(y,t) oscillate at the gyro- 
frequency it is not surprising that a running time average 
increases the rate at which such terms damp out. 

It is our opinion, however, that since the derivation 
from Vlasov equation is of a much more fundamental nature 
than the ad hoc approach the slow decay on the correlation 
time scale more accurately represents the physics of the 
situation. Furthermore it should be noted that as y-»o 
the time to reach asymptoticity grows without limit in both 
approaches. 


Preceding page blank 
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FIGURE CAPTIONS 

Figure 1. Complex k plane showing deformed contours. Original 

contour ran along real axis and below k=k Q . 

Contours were pushed up, C^ around pole and upper 

branch cut, and down, C , around lower branch cut. 

2 

Figure 2. Plot of F(e,x) and eH(e,T) vs. t for a = 3/4 and 
e = 10" 1 and 10 -2 . 
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